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PREFACE 


This  report  is  written  for  programmers  who  create  the  geographical  databases  for 
computer-generated  imagery  in  simulator  displays.  Assignment  of  colors  to  surfaces  should  be 
the  least  of  their  worries,  yet  they  are  often  uncomfortably  aware  that  the  procedures  they  are 
using  have  grown  up  under  the  pressure  of  getting  jobs  done  and  lack  a  consistent  rationale. 
Information  about  the  colors  of  natural  landscapes  can  be  found  in  the  color  science  literature, 
but  programmers  have  no  time  to  search  this  literature.  Physical  data  needed  for  automatic 
computation  of  those  colors  can  also  be  found,  but  no  one  has  collected  them  in  one  place  or 
summarized  their  meaning  for  computer-generated  imagery. 

The  report  emphasizes  relative  luminances  of  natural  surfaces  in  addition  to  their 
chromaticity.  The  color  assignment  for  a  surface  determines  both  the  chromaticity  and  the  range 
of  lightnesses  it  can  have  in  a  daylight  simulation.  Because  we  have  made  use  of  reflectance 
data  extending  beyond  the  visible  range  into  the  near  infrared,  the  tables  in  this  report  will  also 
enable  programmers  to  give  proper  attention  to  relative  lightnesses  in  scenes  intended  for  night- 
vision  simulation. 

This  work  was  conducted  by  the  University  of  Dayton  Research  Institute  under  Contract 
No.  F33615-90-C-0005,  Work  Unit  1123-03-85.  The  authors  wish  to  acknowledge  the  valuable 
assistance  of  the  following  individuals: 

William  C.  Decker  IV,  now  at  the  Electro-Optical  Products  Division  of  ITT  in  Roanoke, 
Virginia,  who  supplied  the  Lotus  program  from  which  some  of  the  reflectance  data  and  two  of 
the  sensitivity  curves  were  taken, 

Alex  Firdman,  who  helped  extract  these  data  from  the  Lotus  program, 

Raymond  L.  Lee,  Jr.,  Oceanography  Department,  U.S.  Naval  Academy,  who  called  our 
attention  to  earlier  literature  on  landscape  colors, 

Elizabeth  Brummer,  who  helped  set  up  the  work  sheets  for  our  computations,  and 

Mrs.  Marge  Keslin,  who  edited,  proofread,  and  typed  the  final  manuscript. 
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COLORS  IN  NATURAL  LANDSCAPES 


INTRODUCTION 


Visual  displays  of  computer-generated  imagery  play  an  increasing  role  in  the 
training  of  military  and  civilian  aircrews.  Stand-alone  flight  simulators  bought  by 
military  units  or  commercial  airlines  commonly  have  the  ability  to  display  daytime  as 
well  as  night  imagery.  Electronic  networks  have  begun  to  link  real-time  simulators  at 
widely  scattered  locations  in  order  that  military  units  of  different  types  may  join  in 
exercises  using  a  common  geographical  database. 

The  colors  which  appear  in  such  a  database  depend  more  on  the  designers’  pref¬ 
erences  and  local  traditions  than  on  color  science.  It  has  been  generally  assumed  that 
the  modeler  should  aim  for  "realism”  in  the  choice  of  colors  and  that  the  best 
simulation  is  the  one  which  most  closely  matches  the  appearance  of  the  real  world.  No 
experimental  evidence  exists  to  support  or  reject  this  assumption.  Before  such  evidence 
can  be  gathered,  display  colors  must  be  brought  under  fairly  precise  control.  Most 
simulator  users  have  not  yet  achieved  good  color  control,  although  simple  methods  of 
doing  so  are  now  readily  available. 

If  the  goal  is  realism,  then  the  choice  of  color  can  be  guided  by  computations 
based  on  physical  data.  The  chromaticity  and  relative  luminance  of  a  surface  are 
completely  determined  by  the  spectral  reflectance  distribution  of  that  surface,  the 
spectral  energy  distribution  of  its  illumination,  and  the  spatial  relations  of  surface  and 
source  of  illumination.  Writers  on  computer  graphics  (Hall,  1988;  Meyer  &  Green¬ 
berg,  1986)  recommend  that  color  selection  be  based  on  these  physical  data.  Although 
they  also  piovide  references  to  articles  containing  such  data,  their  recommendation  has 
not  had  much  impact  on  real-time  computer-generated  imagery,  where  modelers  seek 
"realism"  of  colors  more  often  by  artistic  judgment  than  by  color  science.  Given  the 
known  facts  about  object  colors  picked  from  memory,  it  is  not  surprising  that  computer- 
generated  images  are  more  colorful  than  natural  scenes.  In  memory,  grass  is  greener, 
bricks  are  redder,  and  the  sky  is  bluer,  and  most  objects  are  remembered  as  having 
exaggerated  saturation  and  lightness  (Bartleson,  I960). 

This  article  supplies  reflectance-based  chromaticity  coordinates  and  relative 
luminances  for  surfaces  prominent  in  natural  landscapes.  Not  surprisingly,  the 
chromaticities  derived  from  colorimetric  computations  agree  well  with  the 
chromaticities  reported  by  Hendley  and  Hecht  (1949),  who  determined  landscape  colors 
by  visual  matching  with  Munsell  samples,  and  by  Burton  and  Moorhead  (1987),  who 
studied  digitized  photographs  of  terrain  scenes.  Hendley  and  Hecht  call  attention  to  the 
fact  that  natural  landscape  colors  have  a  limited  hue  gamut.  "Green  plants  fall  in  a 
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yellow-green  region  varying  from  SSO  nm  to  575  nm  in  dominant  wavelength.1  Earths 
and  dried  vegetation  are  yellow  to  orange-red  (576  nm  to  589  nm).  Water,  sky  and 
distant  objects  are  blue  (459  nm  to  486  nm)."  In  autumn,  "vegetation  covers  not 
merely  the  summer  range  of  green  plants,  but  also  that  of  earths,  and  extends  beyond 
in  the  red  to  the  end  of  the  spectrum. "  The  gamut  is  also  limited  with  respect  to  satu¬ 
ration  or  colorimetric  purity.  Except  for  some  autumn  colors,  most  of  the  colors 
studied  by  Hendley  and  Hecht  had  excitation  purity2  less  than  40%,  declining  to  11% 
or  less  at  a  viewing  distance  of  3.5  miles.  With  increased  viewing  distance,  all  colors 
also  shifted  toward  blue.  At  distances  1.34  km  and  beyond,  Burton  and  Morehead 
(1987)  also  found  this  shift  toward  a  hue  "indistinguishable  from  that  of  the  sky"  or, 
occasionally,  slightly  bluer  than  the  sky. 

These  previous  studies  contain  a  few  measures  of  daylight  luminances,  but  the 
lightness  information  is  not  in  a  form  readily  usable  by  modelers.  Since  modelers  of 
computer-generated  imagery  need  information  about  the  lightness  as  well  as  the 
chromaticity  of  scene  elements,  the  present  article  also  contains  relative  luminances 
computed  from  reflectance  data.  The  luminance  of  a  reflective  surface,  relative  to  the 
luminance  of  a  perfect  reflector  in  the  same  illumination,  is  commonly  called  the  "lu¬ 
minance  factor"  (LF)  of  that  surface.  The  tables  presented  here  contain  LF  values  for 
daylight  (photopic)  vision,  unaided  night  (scotopic)  vision,  and  two  types  of  image- 
intensifiers  currently  used  in  night-vision  goggles.  Used  with  the  equations  supplied 
below,  these  tables  of  chromaticity  and  LF  should  be  helpful  to  modelers  who  wish  to 
achieve  greater  realism  in  simulator  scenes  and  to  ensure  similarity  of  color  among 
simulators  operating  at  different  geographical  locations. 


1Dominant  wavelength  is  determined  relative  to  a  reference  white.  In  Hendley  and 
Hecht  (1949),  the  reference  was  Illuminant  A  (tungsten  light),  with  a  color  temperature 
of  about  2800°  K.  When  both  the  reference  white  and  the  color  in  question  have  been 
plotted  in  a  CEE  chromaticity  diagram,  the  straight  line  joining  these  points  may  be 
extended  beyond  the  color  to  the  spectrum  locus  or  to  the  purple  boundary.  For  a  color 
lying  between  white  and  the  spectrum  locus,  the  intersection  with  the  spectrum  locus 
defines  the  color’s  dominant  wavelength  with  respect  to  that  reference  white.  For  a 
color  lying  between  white  and  the  purple  boundary,  the  intersection  with  the  purple 
boundary  defines  the  purest  color  of  this  hue.  Since  no  single  wavelength  can  be 
identified  for  this  hue,  it  is  customary  to  refer  to  such  a  color  by  stating  the  dominant 
wavelength  of  its  complementary  color,  adding  a  negative  sign.  The  complementary 
color  lies  at  the  intersection  of  the  same  straight  line  with  the  spectrum  locus. 

^citation  purity  is  defined  as  the  ratio,  in  the  CIE  1931  chromaticity  diagram,  of 
the  distance  between  a  color  and  the  reference  white  to  the  total  distance  between 
reference  white  and  the  color’s  dominant  wavelength. 
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SOURCES  OF  REFLECTANCE  DATA 


Table  1  is  based  on  a  source  of  reflectance  data  which  is  well  known  to  color 
scientists.  Between  1930  and  1942,  E.  L.  Krinov  (1953)  obtained  spectral  reflectance 
data  on  a  large  number  of  terrain  surfaces  in  several  geographical  regions  of  the  USSR, 
using  various  types  of  laboratory  and  field  spectrographs  available  to  him  during  that 
period.  His  report,  published  in  1946  from  the  Aero  Methods  Laboratory  of  the  USSR 
Academy  of  Sciences,  contains  370  reflectance  distributions  taken  from  about  150 
different  types  of  surfaces.  Krinov  grouped  these  "natural  formations"  into  8  catego¬ 
ries:  forests  and  shrubs,  grass,  mosses  and  lichens,  field  and  garden  crops,  outcrops  and 
soils,  roads,  water  surfaces  and  snow,  and  buildings  and  building  materials.  Most  of 
them  were  viewed  from  the  ground;  some  were  observed  from  the  air  at  an  altitude  of 
about  300  m.  In  certain  cases  the  same  type  of  surface  was  measured  several  times  in 
order  to  study  variations  due  to  season,  moisture,  sun  angle,  and  viewing  angle.  Each 
measurement  was  recorded  on  a  separate  photographic  plate,  together  with  the  reflect¬ 
ance  from  a  standard  plate  in  the  same  illumination.  "The  development  of  the  many 
thousands  of  spectrograms  was... a  vast  undertaking  in  which  an  entire  collective  of  la¬ 
boratory  technicians  and  statisticians  took  part"  (Krinov,  1953,  p.  80). 

The  wavelengths  from  which  recordings  could  be  made  were  limited  by  the 
sensitivity  of  the  photographic  plates.  To  study  the  visible  region  of  the  spectrum, 
Krinov  used  Ilford  panchromatic  plates  with  sensitivity  from  400-650  nm.  At  those 
times  when  he  also  measured  infrared  (IR)  reflectance,  he  used  one  of  several  types  of 
Agfa  infrared  plates  sensitive  to  part  of  the  range  700-1000  nm.  The  range  between 
650  and  700  nm  was  generally  not  studied,  and  in  many  cases  the  infrared  portion  of 
the  reflectance  file  does  not  begin  below  720  nm.  The  columns  for  unaided  vision  in 
Table  1  are  computed  from  Krinov’s  reflectances  in  the  400-650  nm  range.  The 
columns  for  aided  night  vision  are  computed  from  the  entire  set  of  reflectances  provided 
by  Krinov,  with  values  in  the  gap  (between  650  and  700  nm  or  higher)  interpolated  by 
smoothing  to  connect  the  visible  and  infrared  portions  of  the  data. 

Krinov’s  report  is  available  in  English  through  a  translation  prepared  by  E. 
Belkov  and  published  by  the  National  Research  Council  of  Canada.  Maloney  (1986) 
and  Maloney  and  Wandell  (1986)  have  used  his  data  to  show  that  naturally  occurring 
reflectance  distributions  may  be  described  by  a  set  of  no  more  than  four  basis  functions. 

Table  2  is  based  on  data  supplied  by  William  Decker  and  described  in  a  paper 
written  at  the  CECOM  Center  for  Night  Vision  and  Electro-Optics  (Decker,  1989). 
Early  in  the  1970s,  the  Center  undertook  to  develop  methods  for  predicting  the  field 
performance  of  image  intensifiers.  This  project  required  reflectivity  data  in  the  spectral 
region  from  400  to  1200  nm,  and  some  data  of  this  kind  were  collected  during  that 
period.  In  the  mid  1980s,  Decker  and  his  colleagues  at  the  Center  improved  their 
prediction  methods  by  obtaining  additional  reflectance  data  and  measuring  spectral 
attenuation  coefficients  for  four  relative  humidities  (0,  30,  60,  and  90%).  The 
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Table  1.  Colorimetric  Data  Computed  from  Krinov  Reflectances.  Files  are  numbered  and  named  as  in  Krinov  (1953). 
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Table  2.  Colorimetric  Data  Computed  from  Decker  Reflectances.  Files  are  numbered  and 
named  as  in  Decker’s  Lotus  program. 
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reflectances  were  collected  with  a  Perldn-Elmer  Lambda  9  Spectrophotometer;  the 
atmospheric  transmission  data  were  obtained  from  LOWTRAN  VI. 

Decker  has  not  published  these  data;  they  are  provided  to  users  as  part  of  a 
Lotus  program  which  computes  target/background  contrasts  for  several  types  of  NVDs 
under  a  variety  of  viewing  conditions.  Almost  100  reflectance  distributions  are  supplied 
in  this  program,  but  about  half  of  them  describe  the  reflectance  of  military  vehicles  and 
uniforms.  These  distributions  cover  the  range  from  400  to  at  least  1200  nm  (and  in 
most  cases  to  2000  nm)  at  10-nm  intervals.  Data  in  Table  2  for  unaided  vision  have 
been  computed  from  Decker’s  data  in  the  range  400-700  nm.  Data  for  aided  vision 
with  ANVIS  or  CATS  EYE  goggles  have  beat  computed  over  the  range  400-900  nm 
using  the  sensor  response  curves  which  were  also  incorporated  into  the  Lotus  program. 


COMPUTED  CHROMATICITY  OF  SURFACES  IN  DAYLIGHT 

Figure  1  shows  the  location  in  CIE  1976  Uniform  Chromaticity  Space  (UCS) 
for  7  oi  the  Krinov  and  59  of  the  Decker  surfaces  under  daylight  (D65)  illumination. 
The  surfaces  are  grouped  according  to  type.  These  chromaticity  coordinates  are  based 
on  XYZ  tristimulus  values  computed  from  each  spectral  reflectance  distribution 
according  to  the  equation 

T  =  E  R(X)  Et(X),  (1) 

where  T  is  a  tristimulus  value  (X,  Y,  or  Z),  R(X)  is  the  reflectance  of  the  surface  at 
wavelength  lambda,  and  E,(X)  is  the  relative  spectral  energy  distribution  of  D65  daylight 
at  that  wavelength,  weighted  by  the  appropriate  CIE  1964  large-field  color-matching 
function  x,  "y,  or  z".  Values  of  E,(X)  were  taken  from  Table  IV  (3,3,8)  of  Wyszecki  & 
Stiles  (1982,  pp.  774-775),  where  the  tabled  values  have  been  normalized  to  give  Y  = 
100  for  a  perfect  reflector.  Therefore,  the  Y  tristimulus  value  computed  for  any 
surface,  when  divided  by  100,  gives  the  LF  for  that  surface. 

Chromaticity  coordinates  are  computed  from  the  XYZ  tristimulus  values  by  the 
following  equations: 


=  X/(X+  Y  +  Z), 

(2) 

=  Y/(X  +  Y  +  Z), 

(3) 

=  4 X/(X+  15 Y  +  3Z),  and 

(4) 

=  9Y/(X  +  15Y  -I-  3Z), 

(5) 

where  x  and  y  are  the  chromaticity  coordinates  in  the  CIE  1931  chromaticity  diagram 
and  u'  and  v'  are  the  corresponding  coordinates  in  the  CIE  1976  UCS. 
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Figure  1 

Chromadcity  of  Natural  Surfaces.  Computed  from  Krinov  (circles)  and 
Decker  (triangles)  reflectance  data  and  displayed  in  CIE  1976  Uniform 
Chromadcity  Space. 
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The  region  enlarged  in  Figure  1  appears  again  in  Figure  2,  where  it  can  be 
compared  to  the  locus  of  spectral  colors  and  to  the  color  gamut  available  in  typical  CRT 
displays.  Two  of  the  Decker  colors  (US  Coast  Guard  orange  and  blue)  do  not  appear 
in  Figure  1;  they  are  plotted  in  Figure  2.  It  is  clear  from  Figure  2  that  all  the  colors 
in  Figure  1  lie  within  a  region  easily  achievable  by  full-color  electronic  displays  of  any 
type. 


LUMINANCE  FACTORS  FOR  DAY  AND  NIGHT  VISION 

Color  selection  involves  specifying  not  only  the  chromaticity  coordinates  but  also 
the  relative  luminance  for  each  surface.  Variations  in  luminance  due  to  shading, 
texture,  or  orientation  will  be  computed  in  real  time  by  the  image  generator.  The 
database  color  table  provides  a  starting  point  for  these  computations  as  a  set  of  three 
numbers,  the  "RGB  code,"  which  specify  the  voltages  for  the  red,  green,  and  blue 
components  in  the  absence  of  shading  or  texture  effects.  The  relative  luminance  of 
surfaces  in  the  scene  is  fully  as  important  as  their  chromaticity.  Indeed,  it  may  even 
be  more  important,  since  rapid  detection  of  high  spatial  frequency  content  in  a  scene 
depends  principally  upon  luminance  rather  than  color  differences. 

Changing  the  luminance  of  a  scene  component  requires  changing  all  three 
numbers  in  the  RGB  code,  and  the  changes  needed  are  rarely  simple  proportions  for 
each  number.  When  a  modeler  thinks  of  color  habitually  in  RGB  terms,  it  is  hard  to 
give  proper  attention  to  luminance  differences  in  a  scene.  Thinking  in  XYZ  terms 
makes  such  attention  easy;  since  luminance  information  is  carried  by  the  Y  tristimulus 
value,  relative  Y  values  give  relative  luminances  for  any  group  of  scene  components. 
Adjustments  in  Y  values,  accompanied  by  proportional  adjustments  in  the  X  and  Z 
values,  will  always  adjust  luminance  without  changing  chromaticity. 

This  point  can  be  illustrated  by  considering  some  values  in  Table  1,  based  on 
reflectances  from  the  Krinov  data  set.  If  the  luminance  scale  available  in  a  display 
ranges  from  0  to  100  cd/m2,  the  XYZ  values  in  this  table  will  give  realistic 
chromaticities  and  relative  luminances  without  any  adjustment.  A  granite  structure  at 
31.2  cd/m2  will  have  about  half  the  luminance  of  a  patch  of  dry  clay  (65  cd/m2)  and 
about  1.5  times  the  luminance  of  a  red  roof.  Indeed,  if  the  scene  is  restricted  to  the 
components  in  this  table,  the  XYZ  values  can  be  used  without  adjustment  when  the 
display  can  produce  no  output  greater  than  about  70  cd/m2,  since  the  highest  Y-value 
in  the  table  is  65.0  and  the  chromaticity  for  that  value  is  close  to  an  equal-energy  white 
(x  =  y  =  .33). 

But  displays  are  likely  to  differ  in  their  range  of  luminance  output.  In  order  to 
use  the  available  range  most  efficiently,  the  modeler  may  set  the  luminance  of  the 
brightest  surface  in  the  scene  to  the  Y-value  of  the  maximum  luminance  the  display  can 
produce  at  that  chromaticity.  Then  all  other  Y-values  in  the  daylight  scene  may  be 
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Figure  2 

Spectrum  Locus  and  Purple  Boundary  in  CIE  1976  Uniform 
Chromaticity  Space.  Region  enlarged  in  Figure  1  is  shown  here  as  a 
square  containing  the  locus  of  D6S  daylight.  Triangle  inside  the 
spectrum  locus  represents  the  maximum  color  gamut  for  a  typical  color 
CRT.  Chromaticities  for  the  US  Coast  Guard  Blue  and  Orange  colors, 
included  in  Table  2,  could  not  be  shown  in  Figure  1  and  are  indicated 
here  by  crosses. 
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adjusted  according  to  their  luminance  factors,  given  in  the  column  LF,  in  relation  to  the 
Y  and  LF  of  the  brightest  surface.  After  X  and  Z  values  are  adjusted  proportionally 
to  the  new  Y  values,  the  RGB  luminances  for  each  surface  can  be  generated  from  the 
equation 


(LR) 

LG 

=  M'x 

(X) 

Y 

(LBj 

{2) 

(6) 


where  LR,  LG,  and  LB  are  the  RGB  luminances  and  M1  is  the  inverse  of  the 
chromaticity  matrix  M  describing  the  display  primaries, 


M  * 


xR/yR  xG/yG  xBfyB 
1  1  1 
zRtyR  zGfyG  zBfyB 


(7) 


The  x-,  y-,  and  z-chromaticity  coordinates  of  the  primaries  R,  G,  and  B  at  maximum 

voltage  can  be  taken  as  approximations  to  the  values  of  xR,  yR . zB,  even  though, 

in  fact,  the  matrix  M  will  vary  somewhat,  depending  on  the  voltage  producing  the 
required  luminances. 

Luminance  factors  have  also  been  computed  for  night  vision,  and  these  values 
are  included  in  Tables  1  and  2.  To  obtain  these  factors,  the  relative  spectral  energy 
distribution  for  daylight  was  again  used,  on  the  assumption  that  the  distribution  of  night 
illumination  in  the  visible  range  is  not  significantly  different  from  the  daylight  distri¬ 
bution.  However,  for  unaided  night  vision,  the  daylight  distribution  was  weighted  by 
the  function  V'(X),  the  scotopic  sensitivity  function  for  human  vision,  and  normalized 
to  give  Y  -  100  for  a  perfect  reflector.  The  scotopic  LF  values  for  neutral  surfaces 
do  not  differ  much  from  the  daylight  values;  reds  and  yellows  decline  in  brightness 
relative  to  greens,  as  would  be  expected  from  the  Purkinje  effect. 

Tables  1  and  2  provide  additional  information  on  relative  effectiveness  of  these 
surfaces  in  stimulating  ANVIS  and  CATS  EYE  night  vision  goggles.  These  data  can 
be  used  as  luminance  factors  for  simulating  the  appearance  of  a  scene  when  viewed 
through  these  NVDs.  Note  that  daylight  LFs  in  this  table  range  from  0.04  to  0.65; 
ANVIS  LFs  range  from  0. 1  to  0.9. 


13 


CONCLUSION 


If  flight  simulator  displays  are  provided  with  colors  and  relative  luminances 
closely  similar  to  those  in  the  natural  landscape,  will  pilot  performance  improve?  Will 
pilots  find  the  displays  more  acceptable,  or  will  they,  perhaps,  prefer  displays  in  which 
the  colors  are  "richer  than  life?”  Are  there  some  training  purposes  which  might  best 
be  served  by  departing  from  the  colors  of  nature,  at  least  for  certain  classes  of  objects 
represented  in  the  display?  None  of  these  questions  has  been  addressed  in  this  report. 
Indeed,  no  research  relevant  to  these  questions  has  yet  been  performed,  largely  because 
simulator  users  have  only  recently  begun  to  apply  color  science  technology  to  the 
control  of  display  color. 

Nevertheless,  there  is  a  tacit  assumption--yet  to  be  tested  in  controlled 
experiments— that  a  closer  approximation  to  realism  means  an  improvement  in 
simulation.  The  assumption  seems  particularly  appropriate  to  simulations  which  are 
intended  for  use  in  military  mission  rehearsal.  This  report  presents  the  necessary  data 
for  realistic  chromaticities  and  lightnesses.  If  it  encourages  the  simulator  community 
to  give  such  realism  a  real  trial,  the  report  will  have  served  its  purpose. 
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